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■ We study the evolution of the curvature perturbation on the super-horizon scales 

starting from the inflationary epoch until there remains only a single dynamical degree of 
freedom, presureless matter, in the universe. We consider the cosmic inflation driven by 
a multiple number of the inflaton fields, which decay into both radiation and pressureless 
matter components. We present a complete set of the exact background and perturbation 
equations which describe the evolution of the universe throughout its history. By applying 
these equations to the simple but reasonable model of multi-field chaotic inflation, we 
explicitly show that the total curvature perturbation is continuously varying because of 
the non-adiabatic components of the curvature perturbation generated by the multiple 



""^ kiyoung.choi@uam.es 
^ j gong@lorent z . leidenuniv . nl 
^ dj eong@astro . as . utexas . edu 
^Present address 



infiaton fields throughout the whole evolution of the universe. We also provide an useful 
analytic estimation of the total as well as matter and radiation curvature perturbations, 
assuming that matter is completely decoupled from radiation from the beginning. The 
resulting isocurvature perturbation between matter and radiation is at most sub-percent 
level when the masses of the infiaton fields are distributed between 10~^mpi and 10~^mpi. 
We find that this result is robust unless we use non-trivial decay rates, and that thus, 
in general, it is hard to obtain large matter-radiation isocurvature perturbation. Also, 
by using the 6N formalism, we point out that the inflationary calculation, especially 
when involving multiple infiaton fields, is likely to lose the potentially important post- 
inflationary evolution which can modify the resulting curvature perturbation. 
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1 Introduction 



Currently, inflation [Ij is the leading candidate of resolving many cosmological problems such 
as the horizon, flatness and monopole problems, and of explaining the initial conditions of 
the successful hot big bang universe. A certain amount of quasi-exponential expansion of 
the universe during inflation out of small enough region for maintaining causal communication 
makes the whole observable patch of the universe homogeneous and isotropic as well as spatially 
flat. At the same time, the quantum fluctuations of one or more scalar flelds, the inflaton flelds, 
which dominate the energy density of the universe during inflation are stretched and become 
large scale perturbations of spatial curvature. These curvature perturbations are the seed of 
the large scale structure and the temperature anisotropy of the cosmic microwave background 
(CMB) observed today [2]. Recent precise cosmological observations are all consistent with this 
picture and thus indeed strongly support the inflation paradigm [3], H] . 

In spite of the extreme success of the inflation paradigm, the very inflation relevant for our 
own universe is still veiled in mystery. It is because of the lacking of knowledge on particle 
physics: the standard model of particle physics (SM) is one of the most precise theory of 
physics [5], but also is a limited effective theory of what describes all the interactions in nature 
at the very high energy scales relevant for the early universe. For example, the SM does not 
explain the neutrino mass [6] , the baryon asymmetry ^ , dark matter [8j , as well as the inflaton 
flelds. Fortunately, we have several reasonable and attractive extensions of the SM such as 
supersymmetry [H] and extra dimensional theories [TU] , which presumably will be tested at the 
Large Hadron Collider at CERN [11] in the next couple of years. The extended theories of the 
SM commonly postulate the existence of a large number of scalar flelds other than the Higgs 
field which is the sole but yet undiscovered scalar field of the SM. Therefore, it is very natural 
to consider infiation models driven by these multiple number of scalar fields as a plausible 
realization of infiation in this context [T^ . 

The only way to investigate and falsify infiation models is, at best of our current knowledge, 
to study the primordial curvature perturbation. It is observed to be nearly adiabatic and 
Gaussian, with an almost scale invariant spectrum [3l H]. We naturally expect all of these 
properties for single field infiation models. It can be intuitively understood as follows: as the 
only degree of freedom during infiation is the single infiaton field, basically the fiuctuations of the 
infiaton translate into time difference. This gives rise to perfectly adiabatic perturbations. Also, 
suppressed self interaction of the inflaton fleld required by the observed small deviation from 
the scale invariance ensures different Fourier modes of the fleld fluctuations almost perfectly 
independent from each other. This leads to nearly perfect Gaussianity. Thus, any deviation 
from perfect adiabaticity or Gaussianity can be a strong indicator of the existence of more 
than a single inflaton fleld. According to the most recent Wilkinson Microwave Anisotropy 
Probe (WMAP) 5-year data [Ij, the observational constraints are such that the isocurvature 
perturbation between matter and radiation contributes to the total density perturbation on 
the CMB scales at most about 10%, and the non-linear parameter —9 < /nl'^^ < HI 95% 
confldence level, which conflrms Gaussianity of the primordial perturbation at 0.1% accuracy. 

For inflation models driven by a multiple number of flelds, in contrast, there is no such 
unique predictions. Given a fleld trajectory, there are more than one orthogonal directions into 
which the fleld can be 'kicked', then the trajectory can be changed from one to another, i.e. the 
history of the universe can be entirely changed. Thus we need to keep tracking the history until 
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different trajectories coalesce and tliere remains effectively only one degree of freedom: matter 
domination [13]. Tlierefore in principle it is never sufficient to consider only the inflationary 
phase to make any cosmological predictions, as still a large number of the infiaton fields survive 
after inflation and they can subsequently change the history of the universe after the end of 
the inflation. For example, let us consider the curvaton scenario [H]: in that scenario, one 
or more oscillating light scalar fields become dominant component of the energy density of 
the universe in the homogeneous radiation background. The observable quantities such as 
the power spectrum and its spectral index can be completely different from naive inflationary 
predictions, depending on the energy fraction of the curvaton fields at the moment of their 
decay. These curvaton fields themselves can be interpreted as the infiaton fields which have not 
yet decayed until the radiation dominated era, as the individual infiaton field can easily satisfy 
the requirements any curvaton candidate should do [T5] . 

Therefore, on completely general ground, we should consider the entire evolution of the 
universe until matter domination beyond the inflationary epoch when we consider multi-field 
inflation models. The purpose of this paper is to provide a complete set of background and 
perturbed equations which describe the universe starting from multi-field inflation, and to study 
the evolution of the curvature perturbation. Especially, we are interested in the isocurvature 
perturbation between matter and radiation which can be observationally detected in near future 
and thus may be a powerful probe of the early universe. We also discuss the possibility of large 
isocurvature perturbation. 

The outline of this paper is as follows. In Section [21 we present the background and per- 
turbed equations of the components which constitute the universe. They are the infiaton fields, 
and the decay products of the infiaton fields: radiation and matter. We write the equations in 
the spatially flat gauge which is particularly useful for our purpose. We then give the equations 
of the curvature perturbation associated with individual component as well as that of the to- 
tal curvature perturbation. In Section [31 we present the numerical results based on the simple 
chaotic inflation model. We also supply analytic arguments which helps understand the results. 
In Section [H we discuss the importance of the post-inflationary evolution and the possibility of 
generating large isocurvature perturbation. We present our conclusion in Section 5. 

2 Evolution equations 

In this section we present the evolution equations of both background and perturbed quantities 
which constitute the universe. We assume that the universe initially is dominated by a multiple 
number of massive scalar fields, 0j, i = 1,2,3, ■ ■ ■ A^fieids? which play the role of the infiaton 
fields. The effective potential is taken to be 



but nevertheless most of the contents we present in this section should be applicable to more 
general classes of potential. Here 0j decays into radiation and matter with the decay rates F^*^ 



l^ = \/(01,02,---)- 




Especially, we are interested in the separable potential 
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and Tm respectively, which are fixed by underlying physics, e.g. a specific string compactifica- 
tion. We assume for simplicity that the matter product hardly interacts with radiation even 
right after their production. We consider a spatially fiat universe, with the metric including 
the linear scalar perturbations given by 

ds'^ = -(1 + 2A)dt^ + 2aB^idtdx' + a^[{l - 2^)% + 2E^,j]dx'dx^ . (3) 

Note that two of the four scalar perturbations A, B, ijj and E can be eliminated by specific 
gauge choice. In this paper, we pick the spatially fiat gauge where the spatial metric is set 
to be unperturbed, because the perturbation equations take particularly simple forms in this 
gauge as we will see soon. 



2.1 Background equations 

First we present the background equations of motion of the infiaton fields (pi, including the 
decay into radiation and matter, as well as those of the radiation and matter components 
generated by the decay of the infiaton fields. For phenomenological description of the decay, we 

add an additional friction term T^^^(pi = (vi^^ + Fm ^ 4>i to the background equation of motion 

of (pi so that 

4>i+{3H + T^'^)(Pi + V, = 0, (4) 

where Vi = dV/d(pi. It is more convenient to write the equations in terms of the number of 
e-folds using 

dN = Hdt , (5) 

from which we find 

d"^ _tt2 P + P d 

dt^ dm 2ml^dN' ^ ^ 

with mpi = (SttG)^^/^ being the reduced Planck mass. Here, p and p represent the total energy 
density and pressure of the system and are given by 

P =Pj + Pm + ^Pi, (8) 

i 

P=P^+Pm + ^Pi, (9) 
i 

respectively, where the energy density and pressure associated with (pi are 

P^ = \^' + V^ = ^cPf + V,, (10) 

P^ = \^'.-v. = ^<pf-v.. (11) 



3 



with a prime denoting a derivative with respect to A^. Then Eqs. and become 



i ^ ' 



(12) 
(13) 



and 

i 

We can determine the total energy density and pressure by solving the background evolution 
equations of the components. Using Eqs. (EI) and ([7]), Eq. (jll) becomes 

Now we need the evolution equations of the energy densities of radiation and matter pro- 
duced by the decay of the inflaton fields (pi. They can be obtained from the continuity equation 
of a certain component a, which is derived from the conservation of energy-momentum tensor 

T'^^, = ^Tr;. = 0. (16) 

a 

Here a semicolon denotes a covariant derivative and the individual component T^'^ needs not 
be conserved separately, and is given by 

pa = -^H{pa+Pa) +Qa, (17) 

where Qa denotes the rate of energy transfer related to the a component. In the present case, 
the only energy transfer is by the decay the inflation fields with the rate T^^\ For the inflaton 
(pi, radiation and matter, they are given by 

= -r«02 = -if2r»0f , (18) 

i i 

Q».=Er^^^'^' = ^'Er^:Vf , (20) 



respectively. Then, from Eq. ([T7j) the evolution equations are given by 

p; = -4p, + ff5^r«0f , (21) 

i 

p:„ = -3p^ + if5^r«0f , (22) 

i 

for radiation and matter, respectively. 
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Finally, the Hubble parameter H is determined by the Friedmann equation 



= , (23) 



with p being given by Eq. (fT2|) . By taking a derivative with respect to and using the 
continuity equation 

p' = -3(p + p), (24) 
we can write the evolution equation of H as 



Eqs. (|T5]), (in]), (|22]) and ([25]), supplemented by Eqs. ([12]), ([TS]) and dH]), are the equations we 
have to solve to find the evolution of the background quantities p^, pm and H. Note that 
the total background energy density p evolves according to Eq. fl2^ . 

2.2 Perturbation equations 

In this section, we present the perturbation equations to study the evolution of the curvature 
perturbations. As mentioned before we will work in the spatially fiat gauge where the spatial 
components of the metric perturbation are set to be zero. The reason why we use the spatially 
fiat gauge is as follows: from the gauge-ready form of the 00 component of the perturbed 
Einstein equation in a fiat universe, 

(i, + HA)-^[i, + H {a'E - dB)] = , (26) 

we find that the off-diagonal components of the metric perturbation B and E disappear on the 
super-horizon scales since they appear only in the spatial gradient term. Thus, in the spatially 
flat gauge and on the super-horizon scales, the only remaining metric perturbation A is directly 
related to the total density perturbation 

6p = 6p^ + 6pm + ^ Spi (27) 



by eliminating the spatial gradient term and ip from Eq. ( l26l) . 

^ = 4, (28) 
2p 

which makes us free from solving the equation of motion of A: we can determine A once we 
flnd the total background energy density p and its perturbation 6p. 

We derive the perturbed equation of motion of the energy density pa by linearly perturbing 
the continuity equation, Eq. (ITB]) . we can flnd 

6pa + 3H{6pa + 6pa) = QaA + 5Qa , (29) 
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where we take the hmit of our interest, i.e. in the flat gauge and on the super-horizon scales. 
Now we can use this equation to derive the equation of 5 (pi as well as 5p^ and 5pm- Here one 
important point to write the perturbation equations is that, perturbing 0f = —g^^(pifi(pifi would 
incorporate the metric perturbation 6g^^ = 2A. That is, 

50. - = 2i7V: (s<p'i + ^0:) , (30) 

where we have used Eq. (1281) . We can find the perturbed energy density and pressure of (pi by 
substituting Eq. into Eq. ([29]) and using Eq. ([30D, 

5p, =H^[ (^5(P[ + ^(/.'.^ + VMi , (31) 
5p, =H'(P[ (S(P[ + ^0:^ - V,S^, , (32) 



2p' 

respectively. Thus, using these equations and Eq. (jlj), the background equation of 0j with 
additional assumption that F*^*^ is a constant the equation of motion of 6(f)i on the super-horizon 
scales is 

5'^i + 3H5'(l)i + r(')50i - A^i + T^'^A(Pi + 2AVi + ^ l^,,50,- = , (33) 

j 

or in terms of the number of e-folds A^, 



H 2ml^Hy \2p 



P 
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(34) 

which gives the evolution of 5(j)i coupled to the metric fluctuation. Note that with the potential 
given by Eq. ([2]), we have 

(35) 

Also we have to solve the perturbation equations of the radiation and matter energy densities, 
(5p^ and 5pm- We can find the equations by using Eqs. fl29|) and fl30|) as 



5p'^ + 35p^ -hY, ^[11 (^f ^ + 20:^0:) =0 , 



(36) 
(37) 



for 5p^ and 5pm, respectively. 

Finally, we need to supply 5p and 5p which enter Eqs. (!36|) and (!37|) . These are 

completely determined by solving these equations. 5p and 5p are given by 



(38) 
(39) 
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respectively, so that 




Thus, in the spatially fiat gauge on the super-horizon scaled, the evolution equation of the 
total energy density perturbation 5p follows by perturbing the background continuity equation, 
Eq. (Ell), as 

5p' = -3(5p + 5p), (41) 

with 5p+5p given by Eq. fHOj) . Therefore, the perturbation equations to be solved are Eqs. fl34|) . 
f l36l) and fl37|) . Combined with Eq. fj^Tl) and the background equations, we can completely 
determine the perturbed quantities 5 (pi, 6p^ and 6pm as well as 6 p. 



2.3 Curvature perturbations 

With the calculated background and perturbed quantities, we can write the curvature pertur- 
bations associated with the total energy density and a specific component a, which we denote 
by ( and Ca, respectively. The gauge invariant curvature perturbation ( on the uniform density 
hypersurface is defined by [17] 

-C = ^ + H^-l, (42) 
P 

thus in the fiat gauge, ( and (a which is determined by 6pa in the same as ( are written as 

C-^. (43) 
C« = -%. (44) 

Pa 

From Eqs. fH51) and flH|) . it is straightforward that the total curvature perturbation is a 
weighted sum of all the individual curvature perturbations (aS, i.e. 

c = E7C«- (45) 

a 

Now, writing (^s explicitly, we have 

/ = (46) 

C = (47) 

ct>[ m + 0:^p/(2p)] + vMH' 

(3 + r»//7)0f ' ^ ^ 



^Keeping spatial gradient terms, the equation is written as, in the gauge-ready form, 

V2 



5p + 'iH{5p + 5p)-i{p + p) 
where Vu is the perturbed 3-velocity of the fluid. 



= 0. 
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where we have used 



(50) 



(49) 



Thus, Eqs. (ge]), (gTD 



fj48l) . and the total curvature perturbation 



c 



P'^, + P'm + J2i P'i 
Sp, + 6pm + E^ {^'i M + <P^Sp/iM + VM^/H'} 

4p^ + 3p„ + 3//2^,0f 



(51) 



as well as the isocurvature perturbation between any two components [IS] 



<Sa/3 — 3 (Co — Cp) 




) 



(52) 



are completely determined by solving the background equations, Eqs. ( |T5l) . ( 12T|) . ( l22l) and ( l25l) . 
and the perturbation equations, Eqs. (IMI) . (136|) and (137|) . These solutions describe the evolution 
of the curvature perturbations given by Eqs. (l46l) . (H71) . (HH!) and (ISTj) in the universe composed 
of the inflaton fields (pi and their decay products, radiation and matter. Note that although 
(m and Q may meet singularities when the denominators become zero, that of ( remains 
always finite and thus ( is well defined throughout the evolution of the universe. 

3 Evolution of the curvature perturbations 

In this section, we use the equations we derived in the previous section and explicitly study a 
universe initially filled with a multiple number of the inflaton fields until pressureless matter, 
which is produced by the decay of the inflatons, becomes the most dominant component. For 
simplicity, we consider the case of multiple chaotic inflation with the potential 



We first solve the equations given in the previous sections numerically in Section 13. H then esti- 
mate the resulting curvature perturbation along with the isocurvature perturbation analytically 
in Section [3l2l 

3.1 Numerical results 

In this section, we will present the details of the numerical implementation we adopt in order 
to solve the equations we presented in the previous section. 




(53) 
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Figure 1: The probability density functions, p{m), of three different mass distributions: 
Marcenko-Pastur distribution, uniform distribution in hnear scale, and uniform distribution 
in log scale. We use rh = 10~^mpi and j3 = 0.5 for the Marcenko-Pastur distribution , and 
other distributions are uniform between mi and m^^^^^^ given by the Marcenko-Pastur distri- 
bution in linear and log scales, respectively. 
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3.1.1 Parameters 



We consider the case of iVficids = 1000 inflaton fields whose masses are assigned by three different 
distributions: 

(i) Marcenko-Pastur distribution 

The probabihty density function of the Marcenko-Pastur distribution [191 120] is 



= o 0-2 2 ' (54) 



for fn?[l — y/^Y < m? < m^(l + v^)^? ^ind zero otherwise. Here, rh is the average 
of the mass squared, i.e. (m^) = and /5 is a model dependent parameter which 
quantifies the broadness of the distribution. In the calculation, we set rh = 10~^mpi and 
f3 = 0.5, and the inflaton field masses are distributed between mi = 1.69959 x 10~^mpi 
and mjVfl^ids = 2.97069 x 10~^mpi. 

(ii) Uniform distribution in linear scale 

We distribute the mass uniformly in linear scale between mi and ^Affleids °f ^^^^ (i)- That 
is, the probability density function is 

p(m) = (55) 

when rriMfi^^^^ < m < rrii, and otherwise zero. 

(iii) Uniform distribution in log scale 

Finally, we distribute the mass uniformly in log scale between mi and ^Affleids '^^^^ (i)- 
In this case, the probability density function is 

p(logm) = - 1 (56) 

log mi - log m^^^,^^ 

when mTVflpids < m < mi, and otherwise zero. 

We show the probability distribution function for each distribution in Figure [H Note that we 
do not randomly distribute the mass. Instead, we invert the cumulative probability distribution 
function F{x) = p{x')dx' , and assign the mass of z-th field by 

with i = 1,2,3, ■ ■ -A^ficids, so that (pi and (pNti^^^ become the heaviest and the lightest fields, 
respectively. 

Since we do not know the exact decay rate, we use F*^*) proportional to the cubic power of 
the mass of the corresponding inflaton field. In addition, in order for the matter component to 
remain sub-dominant since the big bang nucleosynthesis (BBN) until matter-radiation equality. 
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the decay rate to radiation has to be at least 10 larger than that to matter [2T]. Therefore, 
we choose 



rW =C^ , (58) 

r« =io-^r« , (59) 

where we require C > 10^: this condition comes from that the minimum decay rate we can think 
of is that of gravitational decay, where the decay rate is proportional to the cube of the mass, 
thus r^;^ > mf/mlg. Note that by this difference we are assuming that matter component is 
already non-relativistic when it is produced by the decay of the inflatons. 

Therefore, we can choose the constant C arbitrarily while the highest decay rate, T^^\ 
remains much smaller than the lightest mass, mj^f^^^^^, and we use C = 10^ in this paper: in the 
current case, we require 



r^'^ = Czrf « ^iVH.M. , (60) 

which leads to 



C « ""^"-"-""P' = Cent = 6.05097 X 10« . (61) 

As long as the infiaton fields decay while they are all rapidly oscillating, the final curvature 
perturbation is almost the same irrespective of C, i.e. as long as C < Ccrit- Then, how small 
should C be compare to Ccrit? We have tested that for A^fieids = 100 case with the Marcenko- 
Pastur mass distribution of m = 10~^mpi and P = 0.5. In this case, Ccrit = mioorripj ml ~ 10^°. 
By comparing C = 10'', 10®, 10^ and 10^°, we confirm that the resulting curvature perturbation 
remains almost same in the sub-percent level while C < Ccrit- Especially, when C < 10® ~ 
10~^Ccrit, C is almost exactly the same. Therefore, for A^fieids = 1000, m = 10~^ cases, our 
choice of C = 10^ ~ 10~^Ccrit is safe. The differences basically arises because the inflatons 
decay slowly for lower value of C. We show the results in Figure [2l 

3.1.2 Initial conditions 

In this section, we denote the initial value of a quantity by a subscript '(ini)'. For simplicity, 
we assume that the initial field value is the same for every field, and is set to make the number 
of the e-folds at the end of inflation equal to about 60 with the estimatioE0 



'i^i(ini) 

4mpi 



(62) 



We assume that the slow-roll condition is satisfied initially, and set initial field velocity as 

Vi _ "^i0i(ini) 



I / 



*(ini) ott2 ott2 

"^-"(ini) '^-"(ini) 



(63) 



^This also helps computationally, since too small C requires infeasible amount of computing time as the 
decay of the infiaton fields happens very slowly. 
''See, however, Section H?T1 
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Figure 2: Comparison of density parameter (left) and the total curvature perturbation (right) 
for four different decay rates: C = 10^° (red, sohd line), 10^ (olive, dotted line), 10^ (dark green, 
dashed line), and 10'^ (blue, dot-dashed line), where T^^ = Cmf/mp^. We use 100 infiaton fields 
whose masses follow the Marcenko-Pastur distribution with rh = 10~^mpi and (3 = 0.5. In the 
left panel, thick lines show and thin ones Qif, for each case. As large C increases the decay 
rate, the infiatons decay faster then. In the right panel, all lines show the corresponding total 
curvature perturbations. The resulting curvature perturbations remain almost same in the 
sub-percent level when C < Cent = '"^ioo"^pi/"^i- 
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Here initial Hubble constant, i^(ini), is set to satisfy the Friedmann equation 



tt2 ^ 

'^-"(ini) = Pm(ini) + P7(ini) H ^ 0i(ini) + / . "^i 



i(ini) 



(64) 



As any pre-existing matter or radiation energy density will be exponentially diluted away, we 
set them initially zero, i.e. 

Pm(mi) = P7(mi) = . (65) 

The above equations completely specify the initial conditions for the background evolution. 
Now, we set the initial conditions for the perturbations. 

To compare the numerical results with the analytic estimates, which are performed in the 
spatially flat gauge and in the longitudinal gauge (see Section 13. 2p respectively, we match the 
field fluctuations in these two gauges by using a gauge invariant variable [22] 



Qi = 6(t>i + (p'-ip 



so that 



(66) 



(67) 



where the superscripts '(-F)' and '(L)' denote the spatially flat and the longitudinal gauges, 
respectively. Note that the terms in the square brackets on the right hand side constitute the 
curvature perturbation Q, which makes sense since in the spatially flat gauge simply ip^^^ = 0. 
The initial metric perturbation in the longitudinal gauge ipl^^^, which is denoted as $ in Sec- 

For computational 



^(ini)' 

tion 13. 2^ is obtained from Eqs. (I77j) and (1751) along with Eqs. fl7^ and 
^i(ini) typical value 



convenience, 



^i(im) 



(ini) 



27r 



(68) 



We set the velocity of initial perturbation by imposing the slow-roll condition:under the slow-roll 
approximation, Eq. fl33|) and A = 5p/{2p) are approximated as 



0. 



A 



2p 



respectively. Therefore, we set 

-1 



i(ini) 



(69) 
(70) 



(mi) 



+ 



(ini) 



(nu) 



(2m^0i(ini) + r«i7(i„i)</.'(i^i)) 



(71) 



l + r»/(3if(ini)) 

As we assume that there is no matter or radiation energy density to be perturbed, again simply 

^Pm(ini) = 5p7(ini) = . (72) 
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3.1.3 Evolution 



We solve the system of the coupled ordinary equations by using the eighth order adaptive 
Runge-Kutta method [23]. For numerical efficiency, after the end of inflation, we drop a field 
from the dynamics if the energy density of that field becomes smaller than 10~^ of radiation 
and matter energy densities. That is, (pi is dropped from the dynamical equations when both 

ni<10'^n^, (73) 
fi, <10-^fi„, (74) 

are satisfied. We have also checked different criteria, but the results hardly change. This drop- 
out is valid in our example since heavier fields have smaller energy densities during oscillation 
and decay earlier, so that they do not dominate the energy density later again. 



3.1.4 Results 

Figures [3], H] and [5] show the results of our numerical calculations. In these figures, we show 
the evolution of the density parameters (top panel) and the curvature perturbations (bottom 
panel) of the heaviest (green, dotted line) and the lightest (red, dashed line) infiaton fields, 
the radiation (dark blue, dot-dashed line), and the matter (black, dots-dashed line). The total 
curvature perturbation and the density parameter are also shown with orange solid line. The 
end of inflation, when the slow-roll parameter 

becomes unity, is also shown as a thick vertical line around = 60. In order to analyze 
better the evolution of the curvature perturbation after inflation, we magnify the curvature 
perturbation between 60 < < 100 in the bottom panel. 

After inflation ends, the decay products of the inflatons start dominating the universe. Since 
we set the decay rate to the radiation 10^ larger than that to the pressureless matter, radiation 
dominate ffist and then matter domination comes next. We set the end of the reheating period 
and the beginning of the radiation dominated era when the radiation energy density occupies 
more than half of the total one, i.e. > 0.5. The temperature of the universe at this time, or 
the "reheating temperature" Trh, is calculated from [2] 



Pt = '^3*T^ , (76) 

where g^, represents the effective degrees of freedom of radiation, which is related to the total 
number of particle species. For definiteness, we set g^, = 10^, and calculate the reheating 
temperature for each cases. The results are summarised in Table [1] 

In the bottom panels of the Figures in this section, we find that the curvature perturba- 
tion of 01 (heaviest field) and 0iooo (lightest field) show spikes near the end of inflation. These 
spikes are due to the oscillation of the fields, as 0- becomes very small near the maximum ampli- 
tude of the oscillation. In addition to that, the matter and radiation curvature perturbations, 
denoted by Cm and respectively, change drastically during inflation. We can attribute these 
features to the fact that the curvature perturbation of each component is ill-defined when the 
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Figure 3: Evolution of the density parameter (top) and the super-horizon mode curvature 
perturbation (bottom) of the universe filled with radiation, matter and 1000 inflaton fields, 
whose mass distribution follows the Marcenko-Pastur distribution with rh = 10~^mpi and 
f3 = 0.5. We show the density parameters and the curvature perturbations associated with 
the heaviest and the lightest inflaton fields (green, dotted and red, dashed lines, respectively), 
radiation (dark blue, dot-dashed line) and matter (black, dots-dashed line) as well as the total 
ones (solid, orange line). We also denote the end of the inflation as a thick vertical line. In the 
bottom panel, we highlight the late time evolution of the total curvature perturbation along 
with matter and radiation ones. For multi-field inflation, the curvature perturbation changes 
on the super-horizon scales during and after inflation, and there exists a slight amount of the 
matter-radiation isocurvature perturbation. 
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Figure 4: The same as Figure [3l except that the masses of the inflaton fields are distributed 
according to the uniform distribution in hnear scale between mi = 1.69959 x 10~^mpi and 
miooo = 2.97069 x IQ-^mpi. 



16 



1.0 



0.8 



5 0.6 



o 

CL 



0.4 



0.2 



0.0 



0.0025 



o 0.0020 
o 



0.0015 



0) 



0.0010 



> 0.0005 
0.0000 



total 



10 



total 0\ iradiation 
I 

I 



^■^ /matter 
\ / 

.■"\ 
/ \ 
\ 



I 



\ 



70 80 
e-folding (N) 



90 



70 80 
e-folding (N) 



100 




100 



Figure 5: The same as Figure [3l but the masses of the inflaton fields are distributed according 
to the uniform distribution in log scale between mi = 1.69959 x 10~^mpi and miooo = 2.97069 x 
10-^1. 
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Marcenko-Pastur law 


uniform linear 


uniform log 




60.8433 


60.9048 


60.8669 


A^RH 


68.0806 


68.2138 


68.2071 


TRH/lO^^GeV 


5.91741 


5.38037 


5.07352 



Table 1: Number of e-folds at the end of inflation and at the end of the reheating period, and 
the reheating temperature. 

denominator p'^ occasionally becomes zero. However, total curvature perturbation Ctotai is very 
well defined throughout the whole evolution. 

We find that, for inflation driven by multiple scalar fields, the total curvature perturbation 
constantly changes on the super-horizon scale contrast to the single filed case, even after infla- 
tion. The total curvature perturbation increases during inflation, and even oscillates when the 
inflaton fields are in the oscillatory phase. In the radiation and matter dominated epochs, it 
remains constant following the curvature perturbation of the most dominant component, but it 
changes again during the transition phase around the moment of matter-radiation equality. We 
also find that while multiple inflaton fields can generate the isocurvature perturbation between 
matter and radiation, its amount is very tiny compared with the total curvature perturbation. 
For all of the three mass distributions, the isocurvature perturbation is only about 0.1% of the 
final total curvature perturbation. The results are summarized in Table El 





Marcenko-Pastur law 


uniform linear 


uniform log 


Cm 


numerical 


2.42887 X 10-3 


2.62327 X 10-3 


2.16824 X 10-3 




analytic 


2.47077 X 10-3 


2.66871 X 10-3 


2.21415 X 10-3 


C7 


numerical 


2.42991 X 10-3 


2.62421 X 10-3 


2.16868 X 10-3 




analytic 


2.47247 X 10-3 


2.67120 X 10-3 


2.21625 X 10-3 




numerical 


-3.13389 X 10-6 


-2.83071 X 10-6 


-1.34026 X 10-6 




analytic 


-5.10323 X lO-'* 


-7.46318 X 10-6 


-6.29211 X 10-6 



Table 2: Comparison between the numerical calculations in Section [3?T] and the analytic esti- 
mations in Section 13. 2[ We compare the final curvature perturbation of matter and radiation, 
along with the isocurvature perturbation between them. Analytic estimation provides rea- 
sonably good agreement: it predict the curvature perturbations about 2% accuracy, and the 
isocurvature perturbation within a factor of 0{\). 



3.2 Analytic estimates 

As we have seen in the previous section the curvature perturbation C is not conserved even on 
the super-horizon scales during and after inflation in the presence of a multiple number of the 
inflaton fields. To the best of our knowledge, however, there is no way to describe the evolution 
of C throughout the entire regime of the slow-roll inflation, the phase of the inflaton oscillation 
and the subsequent radiation and matter dominated epochs. The difficulty arises because the 
inflaton fields are not simple barotropic system which satisfies p^p = p<i){P(i)), i.e. in general the 
equation of state of the inflaton field is a non-trivial function of (p and cannot be written as 
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a simple function of p<^ like radiation {p^ = p^/3) or matter (p^ = 0). Nevertheless, we do 
find the evolution of ( during each phase in the literatures. In this section, by combining those 
results together we try to understand the behavior of ( presented in the previous section. 



3.2.1 During slow-roll inflation 



Ref. provides a set of general solutions on the super-horizon scales under the slow-roll 
approximation with a separable potential by sum such as Eq. ([2]). With this potential of the 
form V = Vii^cpi), the general growing solutions of the metric perturbation $ and the field 
fluctuation 6(j)i in the longitudinal gauge and on the super-horizon scales are given by 



$ =eCi + 



3^2 
V 



(77) 
(78) 



where e is given by Eq. ( l75l) . Here, Ci and di are constants determined at the moment of 
horizon crossing given by 



di 



(79) 



(80) 



with the subscript -k denoting that both C\ and di are evaluated at the moment of horizon 
crossing of the mode with momentum A; = |k| and Cij = di — dj. Note that while Ci is 
dimensionless, di and thus Cij have mass dimension —2. We will denote the parts depending 
on Ci as the 'adiabatic' modes, while those on di as the 'non-adiabatic' modes. Below, we 
will denote them by the superscripts '(ad)' and '(nad)', respectively. Accordingly, the energy 
density perturbation 6pi can be written separately as 



5pi 



_ i ^(ad) 



+ V'l 



/.(ad) 



+ 



, (nad) 



+ V'l 



, (nad) 



Spf^ + Spl 



(nad) 



where V- = dVi/dcpi. Then, using 



d_ 

dt 



(ad) 



<|)(ad) 



(82) 



we obtain the adiabatic part of the curvature perturbation as 

^^(ad) XA(ad) 



^(ad) ^ <^(ad) _|_ j^^ 



Pi 



$(ad) _^ 



^3) 



^The 'decaying' adiabatic solutions are written as $ = C2H/a and Sipi/ij) — — 6*2/0, respectively, with C2 
some constant. But as can be read from these solutions, their contributions to Q exactly cancel each other and 
we can completely neglect the decaying solutions [25] . 
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which is indeed constant. From this we can recover the well known "adiabatic" result 



Cl = 7^,, (84) 

with TZc being the comoving curvature perturbation. The non-adiabatic part, the source of the 
time dependence, is from Eqs. (!77j) and (1781) . 

(.i - ^2 + Y ■ y^^) 

With these adiabatic and non-adiabatic parts, the gauge invariant curvature perturbations can 
be written as 

r (nad) 

(. = _ Ci - ^("^'i) - H^^^ = -Ci + Ci"'''^ , (86) 

Pi 

i ^ 

It is very important to notice that both the adiabatic and non-adiabatic solutions of $ and 
5(j)i are contributing to the curvature perturbation C,. More importantly, the evolution of Q 
is purely due to the non-adiabatic contribution. At this point, it would be instructive to 
illustrate the origin of the name 'adiabatic' and 'non-adiabatic' modes: indeed, defining the 
intrinsic isocurvature perturbation as 



Pi 
3 



V ' dt\ V 



we find that the adiabatic modes of $ and do not contribute to S^. Also, the isocurvature 
perturbation between two fields is given by, using Eq. (!52l) . 

where the adiabatic modes are, again, canceled each other. 
3.2.2 During oscillation and afterward 

After inflation ends, the inflaton fields, in fact some of them already, enter the oscillatory phase 
near the minimum of the effective potential. Ignoring the pre-existing matter and radiation, 
the final curvature perturbations of matter and radiation after all the infiatons decay can be 
written in terms of the curvature perturbations of the scalar fields during the oscillatory period 
as [2S] 

cr^ = E ^^cf^'^ = -Ci + E s^'''"'^^-^ , (90) 

i i 

Ci'"^^ = E = + E r.Cf , (91) 
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respectively. The coefficients Sj and rj are given by 

-1 



r(») 



(osc) 



E 



"1 q(osc) 




1 - 



(92) 
(93) 



'71 



We note that ^ • = ^ • Si = 1 by definition. Here, Vtf^'^'' is the density parameter of the field 
04 at the beginning of the oscillatory phase, and p^i is the radiation energy density generated 
from the decay of each 0j. Then the isocurvature perturbation between matter and radiation 
after infiation is, from Eqs. (|52l). fl90|) and (l9T|) . 



(nad)(osc) 



(94) 



We can further simplify the above expressions in the presence of a large number of infiaton 
fields, iVfieids 3> 1, with not too different masses. In this case, at the beginning of the oscillation 
phase (or equivalently the end of inflation) each inflaton field contributes to the energy density 
with a similar order of magnitude, so we can approximate 



(osc) 



fields 



(95) 



Using the fact that during the oscillatory phase the universe is practically matter dominated. 
I.e. a oc 

^2/3 so that we can set the initial time as t^"^"^) =2/ (3iJ(osc))5 then we find [26] 

p^ _ r?i]i°--Vrw [3/(2r(-))]^/3gy^^^) ^ y^^)"' 
P71 ri')fii°''Vr(iM3/(2r(i))]V3i7j/6^ ~ r«'/' ' 

where we have used Eq. (!95l) . Thus, coefficients and are simplified to 



(96) 



J- m 

n 

n 

j=i+i 



EJ- m 
fU) 



3r0-)/ra)5/3 



4 Hill /r^'^)'^' + 3rV^ /ro')'/' 



arvvr 



(97) 
(98) 



In Figure O we present the coefficients Si and used in the numerical calculations of the 
previous section. 

Finally we make an estimation on (^^^^'^'>^°^'^\ estimate the order of magnitude of (^^('^^'^)(°^'^)^ 
we first note that near the end of inflation Vi ^ Vj and V/ ~ V^', as many fields are very close to 
the minimum. Thus, from Eq. ( !85ll . the first term is comparably neglected and we can estimate 

/•(nad)(osc) 

Si as 







(nad) (osc) 



^-"(end) 



fields 



(99) 
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Figure 6: The coefficients Si (black) and (red) used in the previous section, calculated by 
Eqs. f l971) and fl98|) respectively for the Marcenko-Pastur distribution with different (3 values. 
As we have set the fraction of the matter decay rate the same, Tm /T^^ = 10~^, Sj is just a 
constant equal to N^J^^^ = 10^^. Meanwhile, it is clear that the lighter fields contribute more 
to the resulting radiation compared to heavier ones. 



where the subscript '(end)' denotes the end of inflation. Also we have made a wild guess 
Vi ~ V/Nfieids, i-e. each potential energy of the field 0j at the end of the inflation contributes 
equally to the total potential energy. In the above we have defined Pij as the coefficients of the 
stochastic variable dj and it is given by 



A^ficlds 
9 f72 

(end) 



for i ^ i ■ 



(100) 



'1 - M 



fields J 



for 



fields 



Note that for the specific case of the potential given by Eq. (l53l) . it is known |27j that if(end) ~ 
\/2'^Affleids/3 with ^«7Vflcj^^ corresponding to the lightest infiaton mass. Thus, we can write 
Eqs. (ISD and (JH]) simply as 



(101) 
(102) 



respectively. Note that all the coefficients Ci, di and (3ij are directly related to the character of 
the infiaton fluctuations at the moment of horizon crossing, and Sj and rj to that of the decay 
of the infiaton fields. This allows us to connect the initial conditions and/or parameters of the 
model with the observable cosmological predictions with a reasonable accuracy (see Table [2]). 
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Now we can write the 'final' curvature perturbation as well as the isocurvature perturbation 
between matter and radiation. After all of the inflaton fields decay out, there remain only 
radiation and matter components in the universe. Radiation and pressureless matter are perfect 
fluids with the equations of state of the form p = p{p), and the curvature perturbation of each 
species does not change. However the total curvature perturbation does change depending on 
the energy density ratio. They are expressed as 



c 



(fin) 



4p^ + 3pm 



(fin) 



Ap^ + 3p„ 



(103) 



(104) 



We have tested the analytic estimations on C,m^\ C-y^"^^ and Sm^\ given by Eqs. fllOip . f ll02p and 



flMl) respectively, against the full numerical results done in Section [3111 We show the results in 
Table [2l We find that the analytic estimations we have derived in this section predict the final 
curvature perturbations of both matter and radiation within about 2% accuracy. Note that 
the accuracy of the analytic estimations changes with the number of the inflaton fields: it is 
improved as the number of the inflaton fields increases. While the analytic estimation captures 
the smallness of the isocurvature perturbation, it works within a factor of 0{1). 



3.2.3 Power spectra and spectral indices 

Now we can explicitly calculate the power spectra. The Fourier component of the field fluctu- 
ation is obtained by using 

<50,(k) = ^e,(k) , (105) 



/2A;3 

where ej(k) is a Gaussian random variable which satisfies 

(ei(k))=0, (106) 
(e.(k)e*(q))=5.,5(=^)(k-q). (107) 

Then, substituting Eq. (11051) into Eqs. fl7^ and flHU]) . we can write the final curvature and 
isocurvature perturbations in momentum space as 



E 



mlyi, 2^ 4p, + 3p^ VI, 



e,;(k) 



\/2F 
9 



2 V2F 



E 



E( 



e.(k) 



:i08) 



(109) 



Note that pj, which encodes the post-inflationary evolution, has mass dimension —1. Then, we 
define the power spectrum of a quantity Ak as 



Va = 



(110) 
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from which it is clear that the power spectrum of the correlation between and should 
be given by (AkBk). Then, we can simply write the power spectra of the final curvature 
perturbation, the matter-radiation isocurvature perturbation and the correlation between them 
as 



V 

' cur 

V- 

I IS 

V 

' r.oT 



27K 



i 



2 

Vi 1 , 



3K 



PI ia 



. i 



(111) 

(112) 
(113) 



respectively. Note that it is clear for the single field case that we have si = ri = 1 and pi = 0, 
thus from Eq. (11091) 5*^7 = 0, i.e. there exists no matter-radiation isocurvature perturbation. 
This is physically reasonable: as the origin of both matter and radiation is the same, the 
curvature perturbation of the inflaton component C^^ is transferred to both C,ra and ^7- Also 
note that Pcorr can have any sign, depending on whether C*^^'^) and are correlated or 

anti-correlated. 

With the power spectra given by Eqs. fillip and (11121) . we can easily find the corresponding 
spectral indices. As everything is evaluated at the moment of horizon crossing where k = aH, 
we hav^ 

dlogk ^ Hdt . 



Thus, by defining 



n 



d\ogV 



dlogk 

we can write the index of the curvature power spectrum as 



(114) 
(115) 



Ur 



2e 



E 



V,, 



7 -P. 



V 



H 



(116) 



It is worthwhile to note that when the universe is completely dominated by pressureless matter, 
the coefficient pi becomes very simple and so do the spectral indices. In this case, from Eq. (11081) 
we can easily find 



Pi 



2V 



P^ 2 VI' 



:ii7) 
:ii8) 



^Below, we omit the subscript 7k- to avoid too messy notations. 
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1C 

' curv 


^curv 


-p 

' iso 




^iso 


-p 

' corr 




'^corr 


Fig 


.m 


6.06920 X 10-9 


0.958442 


5.32094 X 10- 


-13 


0.977479 


-2.38536 X 10- 


-11 


0.975164 


Fig 


.a 


7.06238 X 10-9 


0.958300 


1.86991 X 10" 




0.976821 


-4.47317 X 10- 


-11 


0.975055 


Fig 


.m 


4.84587 X 10-9 


0.952297 


1.71055 X 10- 


-12 


0.971534 


-4.76410 X 10- 


-11 


0.969105 



Table 3: The values of the spectra and the corresponding indices of the numerical calculations 
using the approximation of Pij in Eq. fllOOl) . 



where we have assumed that both Sj and /3ij are independent on the moment of horizon exit. 
Thus, the spectral index of the curvature power spectrum is approximated as 



2e 



2m 



PI 



V 



X 



2 



ViVl' ( 3 , E.^A V 



n -1 



(119) 



Note that if we neglect pi and substitute Eq. (11171) explicitly with Eq. (1531) . we find after some 
calculations 



ncurv-l = -2e-^^, (120) 

which is consistent with the known results [20l [27] where the post-inflationary evolution is 
neglected. In the same way the spectral index of the isocurvature power spectral is written as, 
once after inflation, 



nisn - 1 



2e + 



2m 



PI 



V 



VI 



(121) 



We can see that for small non-adiabatic contributions compared to the adiabatic one, the 
curvature parts are mainly determined irrespective of the decay properties of the infiatons. 
However the isocurvature perturbation is sensitive to them through the coefficients Sj and r,. 
In Table El we explicitly show the values of the power spectra and the corresponding indices 
of the numerical calculations we presented in the previous section. As can be read from the 
table, Viso contributes about 0.1% of "Pcurv Thus the normalization is mostly determined by 
the adiabatic contribution (see the next section), and we have chosen the parameters to satisfy 
the current observations on the order of magnitude basis. Also note that the indices ncurv and 
niso are different. 



4 Discussions 
4.1 5N formalism 

How is our result of the curvature perturbation presented in the previous section, Eq. fll03p . 
or equivalently Eq. (11081) . related to the one from the 5N formalism, for example presented 
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in Refs. [201 EZ]? With the potential V = Vi(0i), in the slow-roll approximation we can 
calculate the infinitesimal e-folding as 

dN = -m-fJ2^d<P,, (122) 

i « 

so that the number of e-folds obtained during infiation is 

ri(cnd) r<t>i(cnd) Y- 



*(ini) i >' <Pi 

If we specify the case V = ^j^i0i/2, we obtain 

iV = ^i^-^l^!^ = iVo + F, (124) 
4mp[ 4mpj 

where for a large enough number of (pi it is usually assumed that the slow-roll phase is a good 
enough approximation until the end of infiation and that the infiationary phase persists until 
very small value of 0,, so 4>i[end) ~ 0. Thus the integration constant F is usually taken to 
be at most of 0{1) and can be neglected when we estimate A^. However, F does have an 
important implication: by specifying F, we can identify the point where infiation ends, as can 
be read from Eq. fll24p . Otherwise, we lose the information on which trajectory the field has 
been evolving along. The reason is that the end point completely specifies the whole evolution 
and tell us which trajectory we choose, i.e. which universe we live in. That is, omitting F 
is equivalent to start from a minimum around which every direction is the same. This can 
completely jeopardies the evaluation of the curvature perturbation which, as can be read from 
Eq. 

C = 6N = 6No + Y: ' (125) 

y-f «0i(end) d(pj 

depends on the end point, or equivalently, on the trajectory: we may obtain an abruptly large 
or small ( and hence by choosing different end points. 

Now we return to the two-field case, and estimate the power spectrum taking into account 
only the adiabatic component, i.e. 

C = -Cu (126) 



or in the momentum space 



a = , ^ , [0i,.ei(k) + 0ve2(k)] . (127) 



The power spectrum is then 



^curv = + C) ■ (128) 



167r^mp[ 



Meanwhile, taking only Nq from Eq. (11241) . A^'o = {(pl^^ + 02,^)/4mpi so that iVo^i = 0i_j,/2mpp 



thus 



(129) 
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-p 

' curv 


'^curv 


with Pij 


4.88637 X 10-9 


0.952437 


without Pij 


4.90246 X 10-'^ 


0.952492 


fractional difference 


3.29428 X 10-3 


5.75981 X 10"5 



Table 4: The values of the curvature power spectra and the corresponding spectral indices with 
and without the post-inflationary (3ij terms. We have chosen the mass distribution the uniform 
log one. The fractional changes are very small, and the reason is discussed in Section 14.21 

which coincides with the result above. Thus, the adiabatic contribution is associated with the 
leading contribution A^o which does not take a specific trajectory we follow into account. This 
is contained in F, associated with the non-adiabatic contribution. That is, by taking only the 
'adiabatic' component of (, viz. — Ci, it reproduces the power spectrum which results from 
SNq. Omitting the information on the end point encoded in the other piece which completely 
specifies the whole past evolution may completely spoil the result. Fortunately, in this case the 
results with the integration coefficient F is not too different from the one without F, as there is 
only a single global minimum. However if there are more than two equivalent vacua, omitting 
F leads to vastly different consequence. Further, in principle the total number of e-folds up 
to the moment when there remains practically pressureless matter in the universe should be 
written as 

N= Hdt+ Hdt+ Hdt. (130) 

To correctly predict the resulting curvature perturbation (, we must take all these pieces into 
account and perturb them. However, Eq. (11231) only corresponds to the first term and does not 
care about the remaining two, which describe the post-inflationary evolution. Thus it is clear 
by comparing Eqs. (11 111) and (I129p that, we need to include the post-inflationary evolution 
as well as the integration constant F for proper description of the curvature perturbation. In 
Table HI we show the differences between the "pure" inflationary predictions on Vcuw and ricurv 
and those including the contributions of the post-inflationary evolution. Since it is obvious 
from Eqs. (llOOp and (11081) that during matter domination there would be no difference with 
Si being the same for all i, we boosted the difference by a factor Fj = Fi exp [—{i — l)/s] with 
s = —999/ log 10^, which makes an exponential distribution of Fi such that -Fiooo = 10^-Fi. We 
set Fi = 0.001. 

4.2 Large isocurvature perturbation 

Now we turn to the possibility of large isocurvature perturbation between matter and radiation 
Sm-y- It is well known that in the single fleld inflationary model, it is impossible to generate 
5*^7 unless we invoke additional mechanism after inflation such as the curvaton scenario [TU 
HH |26]- In the present context, 5*^7 arises because of the existence of a multiple number of 
the inflaton flelds. Thus observing 5*^7 itself or any of its possible consequences amounts to 
seeing the effects of the multi-fleld inflation. From Eqs. fl^Ul) and (jH]), while the adiabatic 
part of the solution is the same, only non-adiabatic contribution can lead to different curvature 
perturbations associated with matter and radiation. Thus, we can make the coefficients Si and 
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Vi very different, and/or make (^^^^'^^^°^'^^ very different. 

First let us consider (^(^^'^')(°^'^\ From Eq. flHUl) . we can see 

(i. ~ -\/r' ~ -m-^ (131) 

so that Eq. 0991) is written as 

(nad)(osc) ^ y(j__J_]_ (132) 



We can immediately read two consequences. First, the contribution of (^^('^^'^)(°'^'^) qq^^l be either 
positive or negative depending on the mass: for the most (least) massive field, (^^^^'^^^°^'^^ jg 
maximally negative (positive). This can be easily read from the numerical results of the previous 
section, where the curvature perturbation associated with the lightest (heaviest) field is larger 
(smaller) than the total one. Second, as inflation lasts longer and longer, jg driven 

towards zero as H is, obvious from Eq. fl25l) . monotonically decreasing. This leads us to conclude 
that it would not be easy to make very different from each other within a separable 

potential with similar order of masses. 

We can also try to make the coefficients Si and different for each field. As noted in Ref. [26], 
there is hardly any model independent prediction on the matter-isocurvature perturbation and 
the detail depends on the specific decay rates. The decay rate we have used for the numerical 
computations is essentially 



(133) 



which is directly related to the mass of the field. While it would not be completely impossible 
to contrive a non-trivial decay rate which can lead to large matter-radiation isocurvature per- 
turbation, it is not very likely that the fields whose masses are spread over less than a decade 
(see Eq. for example) have completely different decay rates to generate abruptly large 
matter-radiation isocurvature perturbation. Especially, in our numerical computations in the 
previous section we set all the same, and this leads to indeed negligible effects of 5*^7 during 
matter dominated era. Therefore, most probably it would be very hard to have large enough 
matter-radiation isocurvature perturbation in the context of the current model. 

Thus, we need to break the assumptions of the model we have explicitly investigated, i.e. 
the potential includes interaction terms and is no more separable and/or the decay rate is not 
constant so that Si spans a wide range of numbers as (see Figure [6]). Further, it should 
be noticed that even if Piso is negligible, the effects of the post-inflationary evolution can still 
affect ncurv we can observe today: as can be read from Eq. (11191) . a non-trivial set of Sj can 
make the difference of ncurv an observationally detectable one from the estimate based on naive 
inflationary predictions. 

It is worthwhile to note again that while the formulae derived in this paper are general 
and can be applied to any separable potential, the resulting small isocurvature perturbation 
in our specific example is a peculiar property of the investigated model, where the masses are 
distributed in the similar order of magnitudes. If the difference of masses is large enough, such 
as double inflation [28j, then it is possible to obtain large isocurvature perturbation. 
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We have also specifically assumed that the matter component is decoupled from radiation 
already at the stage of generation. This can be possible only for the extremely weakly interacting 
particles, such as gravitino or axino if they are produced dominantly from the decay [29]. For 
the thermally produced dark matter, which is currently the popular candidate, the curvature 
perturbation of matter is the same as that of radiation since they had common values during 
thermal equilibrium. In this case there will be no or negligible isocurvature perturbation. 

5 Conclusions 

In this paper, we have studied the evolution of the curvature perturbation during and after 
multi-field inflation. The infiaton fields (pi decay into radiation 7 and matter m with the decay 
rates F^*^ and Fm respectively, which are assumed to be fixed by underlying physics. We have 
presented the exact set of equations which describe the evolution of the background quantities 
0j, p^, p and H, Eqs. f|T5l) . fl2Tl) . fl22l) . flMl) and fl25l) . respectively, as well as those for the 
perturbed quantities 5 (pi, 6p^, 6pm and 6p, Eqs. f lMl) . fl5B]) . (1571) and fHT]) . Using them, we have 
presented the curvature perturbations associated with each component, Cm and Q which 
correspond to radiation, matter and (pi respectively, Eqs. ( H6|) . (H71) and (HHl) . as well as the 
total curvature perturbation ( and the isocurvature perturbation between two components Sap, 
Eqs. dSH) and (IS2]). 

We have applied these set of equations to a particular simple model of multiple chaotic 
inflation, with the potential given by Eq. (1531) . In Figures [3], H] and [5] we have present the 
numerical results, and have found that ( is continuously varying not only during inflation but 
also afterward. We have presented an analytic argument to account for the evolution of (, which 
is solely due to the existence of the non-adiabatic contributions of the curvature perturbation 
(^^^^^^ given by Eq. flHSjl . (^j^^^'^^ is responsible for the evolution of ( as well as the non-zero matter- 
isocurvature perturbation Sm^\ Eq. (IMI) . We have estimated Cm, C7 and Sm-y analytically, and 
they are in reasonable agreement with the full numerical results with 0{1) factor, as shown in 
Table d 

Throughout this study, we flnd several important points. It is clearly and explicitly shown 
that the total curvature perturbation ( is not fixed but is constantly varying on the super- 
horizon scale, even after the end of inflation. This indicates that any predictions of multi-fleld 
inflation based only on the inflationary epoch may not capture important information about 
the evolution of the universe: e.g. Eqs. (11 111) and (11291) . This also imphes that until full 
radiation domination by the decay of the inflaton flelds, the inflatons themselves can modify 
C a la the curvaton mechanism. Finally, the isocurvature perturbation between matter and 
radiation, which may be detected by near future cosmological observations, becomes large 
enough only when the decay rates are highly non-trivial. This suggests that large matter- 
radiation isocurvature perturbation, if ever detected at all, will be a challenge to our current 
understanding on the decay process of the inflaton flelds. 
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